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ABSTRACT 

We present the results of an experiment to image the interacting binary star 
(3 Lyrae with data from the Navy Prototype Optical Interferometer (NPOI), 
using a differential phase technique to correct for the effects of the instrument 
and atmosphere on the interferometer phases. We take advantage of the fact 
that the visual primary of (3 Lyrae and the visibility calibrator we used are 
both nearly unresolved and nearly centrally symmetric, and consequently have 
interferometric phases near zero. We used this property to detect and correct 
for the effects of the instrument and atmosphere on the phases of f3 Lyrae and 
to obtain differential phases in the channel containing the Ha emission line. 
Combining the Ha-channel phases with information about the line strength, we 
recovered complex visibilities and imaged the Ha emission using standard radio 
interferometry methods. We find that the results from our differential phase 
technique are consistent with those obtained from a more-standard analysis using 
squared visibilities (V^^'s). Our images show the position of the Ha emitting 
regions relative to the continuum photocenter as a function of orbital phase and 
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indicate that the major axis of the orbit is oriented along p.a.= 248.8 ± 1.7°. The 
orbit is smaller than previously predicted, a discrepancy that can be alleviated if 
we assume that the system is at a larger distance from us, or that the contribution 
of the stellar continuum to the Ha channel is larger than estimated. Finally, we 
also detected a differential phase signal in the channels containing He I emission 
lines at 587.6 and 706.5 nm, with orbital behavior different from that of the Ha, 
indicating that it originates from a different part of this interacting system. 

Subject headings: techniques: interferometric — techniques: high angular reso- 
lution — methods: data analysis — binaries: eclipsing — stars: individual {(3 
Lyrae) — stars: imaging 



Introduction 



One of the major challenges in optical intererometry is recovering images, since fringe 
visibility phases are badly corrupted by the atmosphere on short timescales. A few attempts 
have been made to reconstruct im ages from the squared v i sibility (V'^), an unbiased e stimator 



of the visibility amplitude (e.g., [Baldwin, et al. I Il996l . iQuirrenbach et al. Ill994l ). but in 



most cases, images in optical interferometry are produced from V"^ and closure phases (e.g.. 



Monnier et al. 1120071 ). A closure phase is the sum of visibility phases around a triangle of 



interferometer array elements; the atmospheric turbulence effects cancel in the sum, but only 
a fraction 1 — 2/A^ of the phase information from an array of N elements is recovered. 

The challenge in optical interferometric imaging is to recover as much of the phase infor- 
mation as possible. Here we describe the technique we developed to obtain Ha differential 
phases from NPOI observations of the interacting binary star j3 Lyrae. We were able to 
recover complex visibilities of this source, and, for the first time, obtain Ha images using 
standard radio interferometric reconstruction techniques. 



Differential phases, also known as phase referencing (see IQuirrenbach Ill999l . iMonnier 



20031 for a discussion on the subject) is an interesting and powerful technique. Here one 
needs multiwavelength observations, obtained simultaneously, and a priori knowledge about 
the structure of the source in a region of the observed spectrum. From that knowledge, 
one calculates the expected phases in that spectral region, corrects for the atmospheric and 
instrumental effects, and interpolates or extrapolates the corrections t o the wavelengths of 
intere st. One of the first attempts to use this technique was presented by lVakili et al. I (Il997l . 
19981 ). who used Ha and He I 667.8 nm GI2T observations of P Cygni and C, Tauri to detect 
structures in the line-emitting regions of these stars. However, they did not reconstruct 
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images based on their measurements; their conclusions were based only on the analysis of 
the phases and vis ibilities of the stars. Anothe r appli cation of the differential phase technique 
was presented by lAkeson. Swain &: Colavita I (120001 ) , who used the Keck Interferometer to 
detect variations of phase vs. wavelength in data from the binary star t Peg. An ideal 
class of objects to which the differential phase technique can be applied is Be stars, which 
are composed of a B star usually unresolved by interferometric observations (< 1 mas), 
surround ed by a disk of gas with a diameter of a few mas, which is detected in Ha line 
emission (ITycner et al. II2005I ). 



The eclipsing binary star (3 Lyrae (HD 174638, HR 7106, FK5 705) has been ex- 
tensively studied since it was discovered to be an eclipsing system over 220 years ago 
( IGoodricke fc Englefield Ill785r) . Nev ertheless, it remains one of the most baffling binary 
systems known (see iHarmanec 1120021 for a comprehensive review). The current model con- 



sists of a ~3 Mq B6-8II star (star 2 in Harmanec's notation) which has filled its Roche lobe 
and is losing mass at a rate of 2x10^^ year~^ to a ~13 early B star (star 1), which is 
completely embedded in and hidden by a thick accretion disk. This accretion disk simulates 
a pseudophotosphere, which appears to an observer to have a spectrum of an A5III star. The 
orbit is very nearly circular, the orbital period is a little over 12.9 days, and the period is 
increasing by 19 seconds per year. Th e system lies at a distance of 270 ± 39 pc, as measured 
by Hipparcos (jPerryman et al. 1119971). The spectrum of j3 Lyrae contains at least six spec- 
troscopic line components ( IHarmanec 1120021 ). The three major components are: (a) strong 
absorption lines from star 2 used to derive the radial velocity curve having an amplitude 
of 370 km s~^, (b) "shell" absorption components seen in the hydrogen, helium, and many 
metallic lines, believed to arise from accretion disk, and (c) strong lines of hydrogen and 
helium seen in emission which vary with orbital phase and with time. One of the features 
that makes (3 Lyrae unique among contact binaries is the likelihood that the strong emission 
lines come from a bipolar jet that has formed near the point where material s treaming from 



star 2 interacts with the accretion disk around star 1 (IHarmanec et al. Ill996l ). 



This paper is organized in the following way. In Section 2 we present the observations 
and reductions. In Section 3 we describe the corrections applied to the data in order to 
obtain differential phases and describe the results obtained for j3 Lyrae. In Section 4 we 
describe the imaging of the Ha emission of this source. In Section 5 we obtain information 
about the orbit of this system from the displacement of the Ha emission relative to the 
continuum photocenter. In Section 6 we compare the results obtained from the analysis of 
the Ha results with those obtained using a traditional analysis, and in Section 7 we give 
a summary of our results. 
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Observations and Reductions 



2.1. Interferometry 



Our observations were taken witli tlie NPOI (lArmstrong et al. Ill998l ) on five different 
niglits between 2005 May 09 and 2005 May 26. Tlie observations were made witli two 
spectrographs, simultaneously recording fringes in 16 spectral channels in the wavelength 
range 560-860 nm. Between one and three baselines were measured with each spectrograph. 
The maximum baseline lengths ranged from 18.9 to 53.3 m. In Table[T]we present a log of the 
observations, including the number of scans observed each night, the baselines used and the 
range of hour angles covered by the observations. Figured] presents the (m, f )-plane coverage 
achieved by the observations on each night. For reference, this figure also shows the array 
layout with the stations and baselines used. The observations of (3 Lyrae were interleaved 
with observations of a nearby calibrator star 7 Lyrae (B911I, V = 3.24 mag), which has a 
similar spectral type and magnitude to that of /3Lyrae. Each calibrator or source scan had 
a duration of 30 s. 

Throughout this pa per we wil l use o rbital phases of jS Lyrae based on the latest ephemeris 
of primary eclipses from Kreiner ( 2004 1^. We assumed a single orbital phase for each night, 
corresponding to the midpoint of the observations (Table [1]). This assumption will have no 
effect on the final results of this paper, because our observations usually lasted for ~2 hours, 
which corresponds to ~0. 6% of an orbital period. Co mparing the values obtained from this 
ephemeris with those from lHarmanec fc Scholz I (119931 ) we find that the two predictions agree 
at the 2% level, which is sufficient for the purposes of the analysis presented here. As a final 
check, we compared the predicted times of minima with visual measurements from AAVSoB, 
spread over a period of 2 months around our observations. We used only validated observa- 
tions available at the AAVSO site and got an agreement better than 0.25 days between the 
observed and predicted minima, which corresponds to ~2% of the orbital period. 

The ini tial data reductions foll owed the incoherent and coherent integration methods 
described in [Peterson et al. I (j2006al ). In the case of the incoherent integration we start with 
data that were processed to produce squared visibilities (V^) averaged into 1 s intervals. 
These data are then flagged t o eliminate points with fringe tracking and pointing problems 
( iHummel et al. I Il998l . l2003al ). The fiagged data set is then used to bias correct the V'^ 
data, which are then calibrated. The calibrated V"^ values are later used to estimate the 
amplitudes of the complex visibilities. 



http : / / www . as . wsp . kr akow . pl/ephem/ LYR . HTMl 



www.aavso.org 
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For the coherent in tegr ation of complex visibil ities, we use the technique presented by 



Hummel et al. I (l2003bl ) (see iJorgensen et al. 1120071 for a different algorithm). We treat the 



visibilities in each 2 ms frame (the instrumental integration time) as phasors. To determine 
the phase shift from one fame to the next, we compute the power spectrum of the channeled 
visibility as a function of delay and measure the phase difference between the peak of the 
fringe envelope and the nearest fringe peak. We rotate the phasors to correct for this residual 
phase and average the complex visibilities in 200 ms subscans, which are used to determine 
the visibility phases. In Section 3, we describe other corrections applied to these data in 
order to obtain Ha differential phases for (3 Lyrae. 



2.2. Spectroscopy 

To estimate the strength of the Ha emission in the f3 Lyrae system, we have obtained a 
high resolution spectrum using a fiber-fed Echelle spectrograph at the Lowell Observatory's 
John S. Hall telescope. The spectrum was acquired on 2005 May 19, which corresponds to 
the middle of our interferometric run. The spectroscopic data were proces sed using standar d 



routines developed specifically for the instrument used in the observations (iHall et al. Ill994j ) . 
The final reduced spectrum in the Ha region reaches a resolving power of 10,000 and signal- 
to- noise ratio of few hundred (see Fig. [2]). 

We obtain an equivalent width (EW) of —1.5 ± 0.05 nm for the Ha profile shown 
in Figure [2l Most of the uncertainty in the EW measurement is due to the uncertainty 
associated with continuum normalization, which is estimated to be of the order of 3%. 
To obtain the total emission in the Ha line we also need to correct for the underlying 
stellar absorption line that has been filled in by the emission from the circumstellar material. 



Cote & Waters 


Jl987 


) . Using 


(3 Lyrae from 


Dobias & Plavec 



V)o = —0.05 mag for 



component, and therefore we estimate that the EW of the total Ha emission is —2.2 nm. 
Considering that the NPOI Ha channel has a width of 15 ± 1 nm and the total emission line 
width is —2.2 nm, we derive a fractional contribution from the stellar photosphere (ccnt) to 
the total flux in the channel of 0.87 ± 0.01. 
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3. Differential phases 



The chaUenge of recovering phase information from optical interferometric observations 
results from atmospheric and instrumental effects changing the phase calibration on a time 
scale of a few milliseconds. Fringe tracking can reduce these variations but not to the level 
needed for good imaging. Even with fringe tracking, the solution starts with 1) short expo- 
sures that freeze the fringe motion, 2) software that measures and removes these variations 
and 3) a coherent integration of hundreds of milliseconds (or more) needed to increase the 
signal to noise to a level sufficient for the rest of the data processing. The result is mea- 
surements of the real and imaginary parts of the visibility. The uncorrected phase variations 
during the coherent integration consist of a constant part and zero-mean fluctuations. The 
fluctuations reduce the measured visibility amplitude but have no effect on the phase. Thus 
these visibility data can be processed using closure relationships or self-calibration. 

However, with a small number of stations, there are signiflcantly more baselines than 
closure phases, so a signiflcant fraction of phase information is lost. By contrast, the dif- 
ferential phase technique described here provides good phase information for every baseline. 
The idea is to recognize that the observed phase consists of three terms, viz. the desired 
source phase (po, and offsets caused by the atmosphere and the instrument: 



All terms have implicit wavelength dependencies. The differential phase technique con- 
sists of 1) measuring (pobs at several wavelengths, 2) modeling the three right-hand-side terms 
in Equation [Hand 3) deriving parameters of the model by fltting it to the measured (pots- For 
the NPOI, (pinst is stable and can be estimated by observing calibration stars; no parameters 
need to be flt. The atmospheric term, (patm, is given by 



where a is the differential air path between two beams, n is the refractive index of air, k the 
wavenumber (27r/A), d the vacuum path difference and (pc a wavelength independent phase 
offset (Jorgensen et al., 2006). Three quantities need to be flt: a, d, and (pc- Thus with 
at least four spectral channels, some information about the source can be obtained for each 
baseline. 



<Pobs — <Po + <Patm + (pinst ■ 





(2) 



3.1. Instrumental and Atmospheric Corrections 



The Ha emitting circumstellar disks of Be stars constitute one of the best cases for 
which one can calculate instrumental and atmospheric corrections to the visibility phases. 
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This correction is based on observations of calibrator stars and the continuum channels of 
the star of interest. Since the photospheres of these stars usually do not deviate significantly 
from point symmetry, we can model the continuum as having zero intrinsic phase. Even 
in the case of rapidly rotating stars, in which the stellar image may deviate from point 



symmetry, the phase s hould vary smoothly and slowly with wavelength (lYoon et al. 112007 



Peterson et al. 



2006al lb[). The same principle applies to calibration stars. These properties 



allow one to use the information from the continuum channels to calculate the appropriate 
corrections to the phases and interpolate them over the channel containing the H« line. 

A first example of the instrumental and atmospheric effects on the phases of a partially 
resolved star is shown in Figure [31 Each panel of this figure shows a different scan of 7 Lyrae, 
while the different lines in each panel represent individual 200 ms subscans. Notice that the 
lines of the different subscans follow a similar cubic-like shape as a function of wavenumber, 
with some scatter. This cubic curve can be explained as an instrumental phase, due to 
uncompensated amounts of glass along the different beams. This curve does not change 
significantly from scan to scan, or from night to night. A further inspection of this figure 
shows that, besides this cubic component, one can also see a significant variation of the phase 
among different subscans. These variations are due to the differential air path between the 
two beams, and th e shape of th e curve is due to the fact that (n — 1) has a significant 



dependence on /c^ (lOwens 1 119671 ). One can also see in Figure [3] that the scatter is much 
smaller in the left panel, during the first scan of the night, increasing significantly toward 
later times (right), consistent with the worsening of the seeing conditions through the night. 

In order to obtain differential phases of Be stars, one needs to correct the data for the two 
components described above. For the instrumental phases we use the values determined by 
averaging all the 200 ms subscans of calibrator stars observed throughout the night. Since the 
effects of the differential air path and vacuum delay are additive, one can reasonably assume 
that they should average down to zero over the course of a night. Also, since calibration 
stars are expected to have zero phase in the spectral range covered by our observations, by 
averaging the complex visibilities of the calibrators one is left only with instrumental phases, 
which can then be subtracted from the program star visibilities. 

We demonstrate the procedure used to obtain instrumental phases in Figure |H Here 
we show the phases obtained by averaging the visibilities of all subscans of 7 Lyrae observed 
on a night, 720 in total. We can also show that the other component affecting the phases of 
the subscans is due to differenti a l air p ath between the two beams. Following the derivation 



presented by [Armstrong et al. I (119981 ) . we expand Eq. [2] into a power series around fco and 



rewrite Eq. [T] in the following form: 

(j)obs = (j)0 + (f>inst + Co + Cik + C2fc^. (3) 
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Starting from this equation and assuming that (j)o = 0, one can then solve for the coefficients 
Co, Ci and C2 for each subscan, such that the residuals produce the best fit to the atmosphere 
phases {(pobs — 4>inst) ■ We show in Figure HI as black lines, all subscans after subtracting 
only the atmosphere effects from the observed phases. Notice that these lines follow the 
instrumental one, with an r.m.s. of 0.15° in the red and ~0.5° in the blue. 



3.2. Differential Phases of /? Lyrae 

The correction steps applied to the data of the program star (/3 Lyrae) to obtain the 
phase in the Ha channel are presented in Figure O In order to simplify the procedures and 
reduce the number of calculations, we take advantage of the fact that the atmospheric effects 
constitute an additive phase term, so instead of dealing with individual 200 ms subscans we 
can work with the average of all visibilities in a scan. The left panel of Figure [5] shows the 
average phases of nine scans of f3 Lyrae. In the central panel we show the phases of the 
individual scans after the subtraction of the instrumental phase, shown as a dashed line in 
the left panel. The last step of the process is the correction of the atmospheric effects, by 
fitting a quadratic function of the form of Eq. [3] to the continuum channels for each average 
scan and subtracting the resulting function from all the channels. 

In the procedures described above, the continuum phases of P Lyrae were assumed to be 
zero. Nevertheless, one should expect to see phase variations across the spectrum. /3 Lyrae is 



(Wilson 


1974; 


Harmanec 


2002) 



of iMourard et al. I (119921 ). we calculate that this system, observed at maximum separation 
with a 50 m baseline aligned with the separation vector, should show a continuum phase (pent 
at the Ha channel of ~ 14°, with a slope of 1° per channel. This case corresponds to the 
longest baseline in our data, but most of our baselines were 34 m long or less. For a 30 m 
baseline, (pent at Ha is 1.6° with a slope of 0.2° per channel. 

Although the continuum phase component varies smoothly with wavenumber and can be 
modeled by our atmospheric phase subtraction method, arbitrarily setting this component to 
zero can in principle influence our results. The complex visibility of the Ha channel results 
from the sum of the complex visibility of the continuum and Ha emission line. Setting 
the continuum phase to zero will change the phase of the emission line visibility by (pent- 
For baselines of 30 m or smaller this is not an important effect, given the small continuum 
phases. However, as shown above, in the case of our longest baseline this can change the 
Ha emission phase by as much as ~ 14°. We discuss in Section 4 how we tested the effect 
of this assumption on the imaging of this source. 
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Figure [6] shows the differential phases of (3 Lyrae on all baselines observed on 2005 
May 19. Here we see results similar to those from Figure [5|, with the continuum phases 
around zero, and the Ha and He I channels showing a differential phase. The strongest Ha 
phase is found on baselines with the longest extent along the E-W direction, with very little 
differential phase s ignals on baselines AN O-A CO and AEO-ANO. Comparing these results 
with the ones from iMourard et al. I (119921 ) and iHarmanec et al. I (119961) we get a consistent 
picture for this system. Based on GI2T observations, IHarmanec et al. I (119961 ) concluded 
that the orbit of the binary component should be roughly oriented along the E-W direction, 
since their observations used a N-S baseline and did not resolve the two stars. On the other 
hand, their line observations detected smaller \^^'s for Ha and He I indicating that the line 
emitting region is more extended along the N-S direction. This result is c onsistent with the 
spectropolarimetric observations from (IHoffman. Nordsieck &: Foxlll998l ). who found that 
Ha, He I 587.6 nm and ultraviolet (A < 360nm) emission are polarized along p. a. ~ 74°. 
This result indicates that this emission is polarized by a structure perpendicular to this 
direction. Another result that co nfirms this geometry is a radio nebula along p. a. 156 ± 4° 
detected by lUmana et al. I (120001 ) . 



Comparing our results to those from IHarmanec et al. I (119961 ) suggests a possible con- 
tradiction, since we do not see a strong differential phase on the baselines oriented closer 
to the N-S direction, at least when comparing their signal to that measured on baselines 
oriented closer to the E-W direction. However, our observations can be explained in a way 
that agrees with the previous results. First, the weaker phase signal along the N-S direc- 
tion can be expla i ned i n part as an effect of the baseline lengths of the two observations. 
Harmanec et al. I (119961 ) used a baseline of 51 m, while for our AEO-ANO baseline we have 
a longest projected length of only ~ 25 m along the N-S direction. Our maximum baseline 
length along the E-W direction is roughly twice that, which can explain the stronger dif- 
ferential phase signal along this direction. Second, we interpret the differential phases as a 
displacement between the Ha line emitting region and the continuum photocenters. Since 
the binary system is expected to be oriented along the E-W direction, this would explain 
the stronger Ha differential phases on baselines aligned closer to this direction. 

Further support for this interpretation is presented in Figure [TJ where we show the 
differential phases of baseline AW0-E06 on 2005 May 19 and 2005 May 26. These observations 
correspond to orbital phases of 0.24 and 0.78, respectively. Here we can see that the Ha phase 
signal changes from negative to positive between the two dates, which can be understood as 
the center of line emission and continuum photocenter changing orientation relative to each 
other. 



We would also like to point out another interesting result from Figure [6|, the variation of 
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the He I differential phases near k ~ 1.4 and 1.7 /xm"^, in particular on baselines AEO-ACO 
and ANO-ACO. One can see in this figure that the He I channel phases vary significantly from 
scan to scan, with the sign of the phase changing from positive to negative in a few scans. 
When we compare this behavior to that of the Ha phases we see that Ha does not change 
sign on the same night. This result indicates that the He I emitting region has a structure 
different fro m that of the Ha region. This conc lusion is supported by the spectropolarimetric 
results from lHoffman. Nordsieck fc Fox I (119981 ). who found that He I 587.6 nm and 706.5 nm 
have a different polarization properties from those of Ha and He I 667.8 nm. This difference 
is due to the way these lines are excited, with the 587.6 [k = 1.70 /im~^) and 706.5 nm 
{k = 1.42 fim~^) being produced in higher density regions, possibly in the inner accretion 
disk. 



4. Imaging 

The ultimate goal of our differential phase experiment is the recovery of complex visi- 
bilities of this star in order to image the Ha and continuum emission. In order to do this 
we start by combining the final phase information described above with the calibrated \^^'s 
from the incoherent integration method to calculate the corresponding complex visibilities 
for both continuum and line channels. Since these visibilities have been corrected for instru- 
mental and atmospheric effects, the continuum channels can be exported into a FITS file 
and imaged using the usual radio interferometry techniques. However, in order to image the 
Ha emission one needs to take into account the fact that most of the light in this channel 
originates from the continuum. 

The correction for the continuum contribution to the Ha channel is calcu lated using the 



Cent va lue obtained from the spectroscopic measurements. Rewriting Eq. 1 of iTycner et al 



(120061 ) in a more general form, we have the following relation for the continuum subtracted 
Ha complex visibility: 

Vna = (Vobs - CcntVcnt)/i^ - Cent)- (4) 

Vobs is the observed complex visibility in the Ha channel, including the line and continuum 
components. Vent and Vhq are the complex visibilities of the continuum and Ha components, 
and Cent is the fractional contribution from the stellar continuum to the total flux in the Ha 
channel. We use the average of the two channels adjacent to Ha to estimate Vent at the Ha 
channel. 



The resulting visibilities are converted into FITS files (one per night) and imaged using 
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(Ivan Moorsel. Kemball fc Greisen Ill996l ). The resulting Ha images are presented in 



Figure [H These images were created using pixel sizes of 0.5 mas and a hybrid weighting 
scheme between natural and uniform weights, with robustness index 0. Just for reference we 
show in each panel, as a white cross, the position of the continuum photocenter, measured 
on continuum images restored from data observed on the same night. The position of the 
continuum was always in the center of the image, consistent with the fact that the continuum 
phases were set to ~ 0° in our procedures. Notice that in most cases the Ha emission has a 
shape similar to that of the restoring beam, which indicates that the emission is unresolved, 
or just partially resolved, with a size comparable to that of the beam. 

The most important result to be taken from Figure [8] comes from the comparison between 
the peaks of the Ha and continuum emission. We can see in the different panels that these 
two points do not coincide, and that the displacement of the Ha photocenter changes relative 
to that of the continuum as a function of orbital phase. Since the Ha emission originates 
in the disk around the more massive star, while the continuum photocenter is located closer 



to the less massive one (IHarmanec et al. Ill996l ). we are effectively imaging the orbit of this 
system. 

We also tested the effect introduced in the Ha images of assuming that the continuum 
phase on the longest baseline is zero. To solve for this effect, we added 14° to the phases 
of both Ha and adjacent continuum channels of the longest baseline before subtracting the 
continuum contribution to the Ha channel. We did not try to apply this correction to 
other baselines because of the small magnitude of their continuum phases. We followed the 
procedure described above and obtained new images using the uncorrected short baseline 
Ha visibilities and the corrected long baseline Ha visibilities. We do not find a significant 
change relative to the images obtained setting the continuum phase to zero. The position of 
the peaks of emission changed less than OAa, indicating that setting the continuum phases 
to zero was a reasonable approximation in this case. 



5. The Orbit of (3 Lyrae 



Although we do not resolve the stars in this binary, our results can be used to derive some 
properties of the (3 Lyrae system. First, since this is an eclipsing binary, it is straightforward 
to determin e the orientation of the systern on the plane of th e sky. The inclination of the 
orbit is 86° (iLinnell. Hubeny fc Harmanec lll998l : lLinnell Il2000l ) and its projection on the sky 
can be reasonably well approximated by a straight line. In Figure [9] we show the position 



■^http: / /www. aoc.nrao.edu / aips / 
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of the Ha emission relative to the continuum photocenter for the five nights of data used in 
this paper. These positions were obtained by fitting gaussians to the images, and the values 
are given in Table [21 

Figure [9] also shows the orbit of the system, obtained by fitting a line with instrumental 
weighting (1/cr^) to the data points. We find that the orbit is oriented along p.a.= 248.8 ± 
1.7°. Notice that the point from 2005 May 18 has a much higher uncertainty in the N-S 
direction and deviates a little from the line. The uncertainty is caused by the fact that 
station ANO was not used on this night, resulting in poorer resolution along this direction. 
We also point out that the fitted line does not cross the origin, which can be attribute to the 
fact that our observations cover only one half of the orbit, from orbital phases 0.17 to 0.78. 

The reliability of the orbital orientation we obtain for the (3 Lyrae system can be checked 
through the comparison with other measurem ents available in the literature. The radio 
observations presented by lUmana et al. I (120001 ) detected a jet-like nebular structure oriented 
along p.a.= 156.5 ± 4°. This nebula is perpendicular, within uncertainties, to the orbit p. a. 
derived here, agreeing with the expectation that jets are launched perpendicular to their 
disks. Further evidence confirming the orbital orientation of this system comes from the 



compa rison of our results with the spectropolarimetric results of iHoffman. Nordsieck fc Fox 
(119981 ). They predict that th e jets should be oriented along p.a.~ 163.5°, similar to that 



found by lUmana et al. I (120001 ) and roughly perpendicular to the orbital orientation we found. 



Finally, we compare our measurements of the displacement between the Ha and contin- 



Harmanec et al. 


1996; 


Harmanec 


2002: 


Linnell 


2000). 


as that of 


Harmanec et al. 


( 


1996) 


, where the orbit is in 



the most massive star (the one which is currently gaining mass) and the X-axis is defined 
by the line joining the two stars. The analysis of the radial velo cities of this sy stem shows 



that t he eccentricity of the orbit is (IHarmanec fc Scholz Ill993l ). According to iHarmanec 
( 120021 ). the separation between the centers of the two stars is 58.52i?0. For Am = 0.9 mag 
the position of the photocenter is at X = AO.GARfT-, and Y = 0. Based on the analysis of 



the velocity curve of the Ha emission IHarmanec et al 



originates in the disk {X = 4.12i?0 and Y 



(119961 ) found that this emission line 
4.97-R0), in a region offset from the center of 



the system, possibly related to the point of impact of the gas in the disk. 

We used these parameters to calculate the orbits of the Ha and contir iuum photocenters, 

and co nverted them to positions in the sky using the inclination derived by lLinnell. Hubeny fc Harmanec 
( 1l998l ) (i = 86°), the orientation derived by us (p.a.= 248.8°) and a distance of 270 ± 39 pc. 
The resulting Ha orbit relative to the continuum photocenter is shown in Figure [TOl Here 
we can see that the orbit has the correct orientation, but has slightly larger semimajor axis 
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than the one measured by us. The maximum calculated separation between the Ha and 
continuum is 0.623 mas, which is ~ 34% larger than the value we measured (0.466 mas at 
orbital phase 0.24). 

The difference between the observed and predicted orbital semimajor axis can be due 
to several effects. First, the uncertainty in the distance of this system is fairly large (39 pc). 
Increasing the distance of this source to 309 pc (corresponding to one standard deviation 
in the distance), the apparent size of the orbit is reduced by 14%. To illustrate this effect 
we show in Figure dU] (as a dotted line) the orbit obtained using a distance of 309 pc. The 
uncertainty in the Am of this system can also introduce errors in the size of the orbit. If we 
assume Am = O.Smag instead of 0.9 mag, we get that the photocenter is at X = 39.48i?0 
and Y = 0, which reduces the orbit size by ~ 3%. 

Another important source of error in our measurements is the correction of the contin- 
uum contribution to the Ha channel. Factors that contribute to errors in this correction are 
the uncertainty in the wid th of the channe l and t he EW of the Ha absorption component. 



In this paper we used the ICote fc Waters I (119871 ) calibration to correct for the underlying 



Ha absorption. However, their calibration is based on main sequence stars and may be 
overestimating the correction for (3 Lyrae, where the brightest component is a B6-8II star. 
We tested this hypothesis by using the extreme assumption that the EW correction for the 
underlying absorption to Ha is 0, which gives Ccnt = 0.91. Using this value, we obtained new 
images and measured the position of Ha emission relative to the continuum photocenter. 
The new measurements are shown in Figure [TT], where we can see a result opposite to the 
one in Figure [TOl In this case the observed semi major axis is 0.8 mas at orbital phase 0.24, 
~ 28%, or larger than the predicted value. 



6. Comparison Between Differential Phase and V'^ Results 

As a way to check the results obtained with the differential phase technique, we use 
the calibrated continuum V"^ and closure phases to obtain information about the binary 
system. These measurements are used to compute the separation p and position angle 9 of 
the binary for each night. The best fit values, which were obtained for Am(670 nm) = 0.9 
mag and assumed to be the same at other wavelengths, are given in Table [3] and plotted 
in Figure [T^l This figure shows that the maximum separation between the two stars is of 
the order of ~ 0.9 mas, although with large uncertainties. Part of these uncertainties are 
related to the fact that our observations do not fully resolve this system, as can be seen in 
Figure [T3l Here we present, for each night, the ^^'s of one scan observed with the longest 
baseline (AW0-E06, 53 m), and the corresponding best fitting models. We selected the scans 
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within 30 minutes of meridian crossing, which is close to the highest resolution achieved with 
this baseline. This figure shows some of the limitations of our measurements. We observe 
higher \^^'s on the night of May 9, when the system is at orbital phase 0.47, and lower \^^'s 
on the other nights, when the stars are closer to their maximum elongation. However, we do 
not see the characteristic cosine wave signature of a binary, which makes the determination 
of the separation of the two stars uncertain. 

When comparing the results obtained from the analysis with the Ha ones, we should 
keep in mind that the former correspond to the separation between the photocenters of the 
two stars, while the latter correspond to the separation between the Ha and continuum 
photocenters. Consequently, we need to convert the Ha measurements to the same reference 
frame as the continuum measurements. Following the discussion presented in Section 5 we 
have that the two stars are separated by 58.52 Rq, the continuum photocenter is located at 
40.64 Rq from the most massive star, calculated assuming Am=0.9 mag, and for simplicity 
we can say that the Ha photocenter is located at 4.12 R© from this star. These numbers 
indicate that the measured Ha separation relative to the continuum photocenter at orbital 
phase 0.25 corresponds to only 62% of the real separation between the photocenters of the 
two stars. Converting the measurements presented in Figures dU] and [TTl we have that the 
observed Ha orbits correspond to separations in the range 0.75 to 1.28 mas between the 
two stars (the range of values corresponds to Ccnt =0.87 and 0.91, respectively). Considering 
the uncertainties, these values are in very good agreement with those from the analysis 
presented in Figure [T^ 

Another way to confirm the Ha results is by comparing the orientation of the orbits 
derived from the two methods. A least squares fit to the data presented in Figure [121 using 
instrumental weighting, gives a major axis orientation along p.a.= 256.1° ± 17.7°. This value 
is slightly different from the one obtained from Ha measurements, p. a. = 248.8° ± 1.7°, but 
they are still in agreement, given the uncertainties. 



As a last check we used the method described by iTycner et al. I (l2005l . l2006l ) to deter- 
mine the size of the Ha disk. We used the calibrated Ha V^'s, corrected to eliminated the 
continuum contribut ion from the binary star to this channel. Following the results from 



Tycner et al. I (120061 ) we decided to fit the Ha V s with a gaussian model, which results in a 



disk with HWHM=0.60±0.10 mas (Figure [T^. The gaussia n HWHM is comparable to the 



Roche lobe radius of ~0.52 mas (30.3 R© iHarmanec 1120021 ). and consistent with the orbit 



obtained using the two methods described above. 
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7. Summciry 

In this paper wc presented the development of a differential phase technique using NPOI 
data. We described the methods used to correct the instrumental and atmospheric effects 
on the phases of /? Lyrae, which allowed us to recover the complex visibility of the Ha 
channel. We also described the procedure used to correct the continuum contribution to 
the Ha channel and successfully imaged the line emission of this system using standard 
interferometric techniques. We found that the major axis is oriented along p. a. 248.8°, 
consistent with previous spectropolarimetric measurements, radio and optical interferometry 
results. From the comparison between the observed position of Ha relative to the continuum 
photocenter and values obtained from models, we find that our measurements indicate a 
semi major axis smaller than that predicted by the models. We suggest that the most likely 
causes for this discrepancy are a wrong distance, brightness ratio between the 2 stars and 
the uncertainty in the correction of the continuum contribution to the Ha channel. We also 
found that the results obtained using the technique developed here are consistent with those 
obtained from a more traditional analysis that uses only measurements. 

The technique presented here represents a major step in the process of obtaining images 
from optical interferometric observations. Due to several limitation, the most important 
being the rapid time scales in which the atmosphere changes, the phases of a source wrap 
around too fast. Consequently, this information is lost and one cannot apply the usual image 
reconstruction techniques employed in radio interferometry. By making some assumptions 
about the structure of the source, our technique allows us to recover this information and 
obtain images for the emission line channels. In the case of j3 Lyrae this allowed us to measure 
the orientation and size of the orbit using the relative location of the Ha and continuum 
photocenters. We also found that the Hel emission has a different structure from that of 
Ha. 

Besides f3 Lyrae we are also applying this technique to other Be systems. This will allow 
us to image their circumstellar disk, map structures and study in better detail the accretion 
flow and effects such as the ionization of the disk by a binary companion. 

The work done with the NPOI was performed through a collaboration between the 
Naval Research Lab and the US Naval Observatory, in association with Lowell Observatory, 
and was funded by the Office of Naval Research and the Oceanographer of the Navy. We 
thank the NPOI staff for the careful observations that contributed to this work. C. T. thanks 
Lowell Observatory for the generous time allocation on the John S. Hall 1.1 m telescope and 
thanks Wes Lockwood and Jeffrey Hall for supporting the Be star project on the Solar-Stellar 
Spectrograph. This research has made use of the SIMBAD literature database, operated at 
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CDS, Strasbourg, Prance. We acknowledge with thanks the variable star observations from 
the AAVSO International Database contributed by observers worldwide and used in this 
research. We would like to thank Christian Hummel for the availability of the data reduction 
package OYSTER. 
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Fig. 1. — The u — v coverage of the different nights. The black points show the continuum 
channels, while the green ones show Ha. The observing date is presented in the top right 
corner of each panel. The bottom right corner of the figure shows the array layout. The 
astrometric stations are shown as black squares and E06 is shown as a black dot. We also 
show the baselines used all nights as solid lines and the baselines that were used only on 
some nights, those involving AN, as dashed lines. The maximum baseline lengths can be 
found in Table [H 
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Fig. 2. — Normalized high resolution Ha profile of (3 Lyrae. 
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Fig. 3. — Coherently integrated pliases of 3 scans of 7 Lyrae, one per panel, observed on 
2005May09 with basehne AW0-E06. Each hne corresponds to a different 200 ms subscan. 
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Fig. 4. — Comparison between the phases of individual subscans from Figure 3, corrected for 
atmospheric and vacuum delay effects, and the instrumental phase (white line), obtained by 
averaging all subscans of the calibrator star observed in this particular night. The uncertainty 
in the instrumental phase is 0.15° in the red part of the spectrum, growing to ~ 0.5° in the 
blue. 
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Fig. 5. — The different correction steps applied to tlie data of (3 Lyrae. Tlie first panel 
(left) shows the average phases of each scan (solid lines) as for baseline E06-AC0, obtained 
on 2005Mayl9. The dashed line shows the instrumental phases and the dotted vertical line 
shows the location of the Ha channel. The middle panel shows the phases obtained after 
subtracting the instrumental contribution. The right panel shows the differential phases 
obtained after subtracting a quadratic function, which was fitted to the continuum points 
only. The uncertainty in the phases is of the order of 0.3° in the red, 0.5° at Ha, growing to 
1° in the blue. 
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Fig. 6. — Residual phases of (3 Lyrae after correcting for the instrumental, atmospheric 
and vacuum differential delay effects. Each panel shows a different baseline, observed on 
2005Mayl9. The individual lines represent individual scan averages. The vertical dotted 
line indicates the wavenumber of Ha, while the dashed ones show the wavenumbers of Hel. 
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Fig. 7. — Comparison of the differential phases of the longest baseline AW0-E06 (53 m) 
on the nights of 2005Mayl9 (left) and 2005May26 (right). The two observing nights are 
separated by approximately half an orbital period. 
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Fig. 8. — The following 5 figures show the Ha images of (3 Lyrae. Each figure corresponds to 
a different night and are organized in order of increasing orbital phase, indicated at the top 
of the panel. The cross shows the photocenter of the system, measured on the continuum 
images. The lowest contour corresponds to the 3a level. The reconstructing beam size and 
orientation is shown in the bottom left corner of each panel. 
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Fig. 9. — Position of the peak of Ha emission relative to the continuum photocenter. For 
each point we give the position uncertainty, obtained by fitting a gaussian to the image 
of the source. Besides each point we also indicate the orbital phase and the date of the 
observations. The dotted hne was obtained by fitting a least squares line to the data. Since 
(3 Lyrae is an eclipsing binary, this line is a good approximation for the orientation of the 
orbit in the plane of the sky (PA=248.8±1.7°). 
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Fig. 10. — Same as Figure 9, but displaying the orbit calculated using the parameters 
described in Section 5 (solid hne). The dotted line shows the orbit obtained if we assume 
that the system is at a distance of 309 pc, corresponding to a 1 cr deviation in the distance. 
The lines connecting the points to the orbit indicate the expected position of the star in the 
model. 
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Fig. 11. — Same as Figure 10, but displaying the Ha positions measured on images created 
using a larger correction for the continuum contribution to the Ha channel {ccnt — 0.91). 
For simplicity we show only the model where the binary is at a distance of 270 pc. 
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Fig. 12. — Separation of the two binary components, obtained by fitting the continuum 
data. For each point we give the uncertainty in the position of the separation of the 
two sources. Besides each point we also indicate the orbital phase and the date of the 
observations. The dotted line was obtained by fitting a least squares line to the data, which 
gives the orientation of the orbit in the plane of the sky (PA=256.1±17.7°). 
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Fig. 13. — Calibrated V^'s of /3 Lyrae observed with the longest basehne, AW0-E06 (53 m). 
The respective nights and orbital phases are indicated inside the panels. In order to make 
the visualization of the plots easier each panel shows only one scan, corresponding to the one 
observed within 30 minutes of local meridian crossing. The solid lines show the best fitting 
model for the corresponding night. 
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Fig. 14. — Calibrated Ha visibilities from the Ha channel. The dashed line shows the 
continuum after the removal of the binary signature. The dotted line shows the best 
elliptical Gaussian model fitted to the data. 
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Table 1. Observations Log 



Date 


UT Range 


Number of Scans 


Baselines 


Hour A 


ngle 


Orbital Phase 


2005May09 


9.63 . . . 11.90 


5 


E6-AC, AW-E6, AW-AC, AE-AC, AN-AC, AE-AN 


-1.50 . . 


0.79 


0.47 


2005Mayl8 


9.85 . . . 11.80 


10 


E6-AC, AW-E6, AW-AC, AE-AC 


-0.67 . . 


1.28 


0.17 


2005Mayl9 


9.58 . . . 11.75 


9 


E6-AC, AW-E6, AW-AC, AE-AC, AN-AC, AE-AN 


-0.89 . . 


1.29 


0.24 


2005May24 


9.08 . . . 11.63 


8 


E6-AC, AW-E6, AW-AC, AE-AC, AN-AC, AE-AN 


-1.10 . . 


1.48 


0.63 


2005May26 


8.98 . . . 11.36 


10 


E6-AC, AW-E6, AW-AC, AE-AC, AN-AC, AE-AN 


-1.04 . . 


1.34 


0.78 



Note. — The maximum lengths of the baselines shown in Column 3 are: E6-AC 34.3 m, AW-E6 53.3 m, AW-AC 22.2 m, AE-AC 18.9 m, 
AN-AC 22.9 m, AE-AN 34.9 m. We always used AC as the reference station. 



Table 2. Displacement of the Ha and Continuum Photocenters 



UT Date 


Orbital Phase 


RA(H«-Continuum) 
(mas) 


DEC{Ha-Continuum) 
(mas) 


2005May09 


0.47 


-0.13 ± 0.03 


-0.02 ± 0.05 


2005Mayl8 


0.17 


-0.44 ± 0.05 


-0.05 ± 0.16 


2005Mayl9 


0.24 


-0.44 ± 0.03 


-0.16 ± 0.04 


2005May24 


0.63 


0.29 ± 0.03 


0.15 ± 0.05 


2005May26 


0.78 


0.39 ± 0.02 


0.18 ± 0.03 



N ote. — The orbi tal phases were calculated based on primary eclipse ephemeris 
from lKreinerl ll2004h 



Table 3. Binary Separations 



UT Date 


Orbital Phase 


P 


e 












(mas) 


(deg) 


(mas) 


(mas) 


(deg) 


2005May09 


0.47 


0.47 


248.1 


0.52 


0.23 


5.0 


2005Mayl8 


0.17 


0.85 


261.5 


1.09 


0.19 


177.9 


2005Mayl9 


0.24 


0.90 


260.0 


0.52 


0.21 


2.7 


2005May24 


0.63 


0.78 


76.7 


0.51 


0.23 


3.8 


2005May26 


0.78 


0.91 


72.5 


0.52 


0.21 


1.9 


Note. — 


Column 1: date of the observations; 


Column 


2: orbital phase; 



Columns 3 and 4; binary separation and position angle; Columns 5, 6 and 7: 
position uncertainty ellipse major axis, minor axis and position angle. 



